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Abstract 



CD 'Formation of anomalous isotope abundances in the atmospheres of chemically peculiar (CP) stars can be explained by 
light-induced drift (LID). This effect is additional to the radiative acceleration and appears due to systematic asymmetry 



Oh 



of radiative flux in partly overlapping isotopic spectral line profiles. LID causes levitation of an isotope with a red-shifted 
spectral line and sinking of an isotope with a blue-shifted line, generating thus diffusive separation of isotopes. We have 



studied diffusion of mercury as a typical well-studied isotope-rich heavy metal. Our model computations show that 
m mercury-rich quiescent atmospheres of CP stars LID causes levitation of the heavier mercury isotopes and sinking 
of the lighter ones. Precise quantitative modelling of the process of isotope separation demands very high-resolution 
computations and the high-precision input data, including data on hyperfine and isotopic splitting of spectral lines, 
i - 1 adequate line profiles and impact cross-sections. Presence of microturbulence and weak stellar winds can essentially 
reduce the effect of radiative-driven diffusion. 
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Q .1. Introduction 

The problem of specifying the physical processes which 
, 1 cause the observed drastic overabundance of several heavy 
metals and anomalous abundance ratios of their isotopes 
t-H dn the atmospheres of chemically peculiar (CP) stars, espe- 
^ cially of chemically peculiar mercury-manganese (HgMn) 
stars, has been on agenda already for several decades. At 
jthe present time almost nobody doubts, that atomic diffu- 
sion is responsible for a large part of the abundance vari- 
^ ations observed in CP stars. About 40 ye ars ago M i chaud 
f- ^ demonstrated in his pioneering papers (iMichaudl . Il970l : 

o 



iMichaud et al. . 19761 ) that for spectral line-rich ions of 
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trace elements the gradient of the radiative pressure, i.e. 
the upwards directed radiative acceleration, essentially ex- 
ceeds the downwards directed gravity. This circumstance 
produces the phenomenon of accumulation of spectral line- 
rich heavy metals in the stellar atmosphere. Further nu- 
merous papers have confirmed the result. 

Stellar evolution models with radiative-driven diffusion 
have been con structed by Michaud and his colleagues fo r 



AmFm stars (IRicher et all. l2000t M ichaud et all 120051 ). 



Popu lation II stars llRichard et all l2002bllal: VandenBerg 
et al. . 2002M and horizontal branch (TI B) stars ( Michaud 
et al.. 120071 . l200£j ). HB stars are very closely related to 
HgMn stars. These models do not include detailed at- 
mospheric modelling and assume chemically homogeneous 
mixed outer region. Radiative accelerations have been 



foun d using so-ca lled "diffusion approximation" for photon 
flux l|Milnd . ll927l ). This approximation is valid in optically 
thick medium of stellar interi ors where r ^> 1. Evq l ution- 
ary computations show (see Michaud and Richeil 120081 . 
Fig. 1) the high trend of cumulation of heavy metals in 
the outer layers of stars. Line-rich metals are swept up 
throughout the stellar interior by upwards directed radia- 
tive push. Radiative acceleration on the extremely line- 
rich low ionization stage ions of heavy metals can exceed 
downwards directed gravity up to 3 dex. 

Stellar evolution models have been relatively successful 
at explaining peculiarities of AmFm stars. However, not 
all abundance anomalies of HgMn and HB stars can be 
reproduced. Isotopic anoma lies observed in HgMn stars 
also remained unexplained. Michaud and Richer (2008) 
suggested that additional separation occurs in the atmo- 
spheric regions. 

Observational evidences of abnormal isotope ratios in 
CP stellar atmospheres have been cumulated during last 
decade owing to the high-resolution and high signal-to- 
noise exposures of the HST and ground-based telescopes 
with similar capabilities. Although the ove rall picture of 
isotope variations is complex (see for review ICowlev et al 



2008), there is a general regularity for heavy elements to 



show overabundance of heavier isotopes and for light ele- 
ments - overabundance of lighter isotopes. For example, 
isotopic abundances of mercury range from the terrestrial 
mixture to the virtually pure 2 04 Hg while overall overabun- 



dance of mercury is 4-5 dex (jDolk et all 120031 ). In spite 
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of numerous attempts, strong isotopic anomalies observed 
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in the atmospheres of chemically peculiar stars have not 
found acceptable explanation yet. Separation of isotopes 
requires some additional physical mechanism, because the 
radiative push is almost the same for all isotopes of the 
particular element. 

A possibility for br eakthrough was elucidate d bv Atu- 
tov and Shalagin (j 19881 ) . who, based on the results of lab- 
oratory laser experiments, proposed a new physical mech- 
anism, called light-induced drift (LID), for explaining the 
isotope separation in the quiescent atmospheres of CP 
stars. 



2. Concept of LID and its main formulae 



The proposal by Atutov and Sha agin dl988f ) inspired 
us ( Aret and Sapar . 20021 ; ISapar et al. . l2008al bfTto investi- 
gate the possibilities to unveil the enigmatic appearance of 
isotope anomalies in the CP stellar atmospheres using LID 
phenomenon. This new mechanism is fully based on the 
asymmetry of spectral lines. It works effectively if spectral 
line profiles of a studied element are systematically asym- 
metrical. This is the case for partly overlapping mutually 
shifted isotopic lines. 

The absorption of the radiative flux, which is asymmet- 
rically distributed in the blue and red wings of a spectral 
line, generates anisotropy of excitation rate of the atomic 
particles participating in thermal motion. For example, 
larger radiative flux in the red wing of a spectral line gives 
more excited particles among those moving downward in 
the atmosphere than among those moving upward. Par- 
ticles in excited states have larger collision cross-sections 
than particles in the ground or lower excitation states and 
thus the mobility of particles in the higher excited states 
is smaller. If particles survive in the excited states until 
the next collision, then their free paths are shorter than 
free paths of the unexcited particles. This means that 
radiative flux asymmetry in the spectral lines triggers a 
diffusion phenomenon. 

The first task was to find the formulae describing this 
phenomenon. We succeeded to reduce the LID to the 
equivalent ac celeration to be added to the usual radiati ve 
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acceleration ([Aret and Saparl . 120021 ; ISapar et al. 

As well known, the usual radiative acceleration due to 
electron transition from lower level I to upper level u of 
ion j can be found as 



^rad 



/ Xj ia u iV(u v , a)T v dv 

rrijC J 



(1) 



Here rrij is the mass of the light-absorbing ion, Xjj is state 
I population fraction, ~kT v is total monochromatic flux and 
a\i is the absorption cross-section in transition I — > u: 



n e 2 f u i 
m e cAi>D ' 



(2) 



where f u i is oscillator strength and Avd is Doppler width 
of the spectral line. Normalized frequency distribution in 



a spectral line is the Voigt function, being the convolution 
of the Lorentz and Doppler profiles, i.e. 



V{u v ,a) 



3 / 2 J (tt„ — y) 2 + a' 



■dy . 



(3) 



The argument of the Voigt function is the dimensionless 
frequency u v = Ai>/ Avf>, its parameter a = T u i/(AttAi/d) 
is the ratio of characteristic widths of Lorentz and Doppler 
profiles. Integration is carried out over the dimension- 
less velocity y = v/vt, where the thermal velocity vt = 
y/2kT/nij. 

LID appears due to difference of collision cross-sections 
in different quantum states. LID efficiency depends on the 
difference between collision rates in upper and lower state 
and also on the probability of particle to stay in the excited 
state until the next collision. 

Based on these conceptual standpoints we have derived 
the following expression for the equivalent acceleration due 
to LID in the transition I — + u of ion j: 



, lid _„ K I Y n dV{u v ,a) 

eq / -Ajja ul J- V av . 
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du v 



(4) 



Coefficient q is the ratio of the mean momentum of the 
atomic particles and the photon momentum: 



rrijVTC 
2hv 



mjVT hv 
2 ' ~c 



(5) 



Due to the large value of q ~ 10, 000, LID can be important 
even in the case of moderate asymmetry of radiative flux 
in the spectral line profiles. Under favourable conditions 
LID may be several orders of magnitude larger than usual 
radiative acceleration in line. Efficiency of LID e is given 

by 

C u A u + C u A u + C u 

where A u is the total rate of radiative transitions from the 
upper state. The first term in Eq.© is the relative dif- 
ference between collision rates in upper and lower states. 
The second term expresses the probability that during the 
free flight after photon absorption de-excitation does not 
happen. Unlike the amplification coefficient q, the value of 
£ can be estimated only with rather low precision. Char- 
acteristic values of e are from about 10~ 3 in the outer 
atmospheric layers of the CP stars to about 1/4 in the 
deep layers. 

Quantum state life-times have been estimated similarly 
to the Kurucz spectrum synthesis program SYNTHE (Ku- 
rucz, Il993h acc ordin g to the formulae given in the book by 



Kurucz et al.l (119741). Radiat ive and Stark damping con- 



stants bv lProffitt et all ( |l999h were used for the lines given 
in their paper. Cross-sections of mercury collisions with 
buffer gas (H + He) particles in the upper state u are larger 
than in the lower state I. These collisions are predominan- 
tely elastic. The collision cross-sections were calculated 
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adopting the quasi-hydrogenic approximation. For ion- 
ion impacts also the additional Coulomb repulsive cross- 
sections were taken into account . Diffusion coe f ficient s 
were calculated according to the IGonzalez et aL ( 1995 ). 
These approximations allow to study general trends of iso- 
tope separation due to LID. More detailed computations of 
cross-sections and quantum state life-times are necessary 
in the future studies, but general picture should remain 
the same. 

The total effective acceleration, acting on atomic par- 
ticle due to both the radiative force and the LID is given 
by a* = aj ad + a^ m . Thus, the effect of the light-induced 
drift can be incorporated by replacing the Voigt function 
V(u v , a) in the expression of radiative acceleration Eq.([T]) 

by 

dV(u v , a) 



W(u v , a) — V(u v ,a) + eq- 



du h 



(7) 



Once again, we would like to emphasize, that the LID can 
be treated as acceleration generated by additional specific 
force. 

Light-induced drift is most effective in the intermediate 
layers of stellar atmospheres. In high layers of atmosphere 
or above it LID efficiency is small because of long free paths 
of the particles, i.e. probability of particle to return to the 
lower state before the next collision is larger. In stellar 
interiors LID is inessential compared to the usual radiative 
acceleration due to other effects: first, isotopic splitting of 
lines of high ions (e.g. Hg III, IV etc.) is so large that 
isotopic lines do not overlap and consequently do not give 
a systematic asymmetry of the lines generating large LID; 
second, radiative flux becomes local and asymmetry of flux 
diminishes. 

3. On the isotopic and hyperfine splitting of spec- 
tral lines 

The efficiency of the LID phenomenon is governed by 
the finest details of spectral line splitting, including their 
hyperfine and isotopic structure, which determines the reg- 
ularities in the spectral line overlapping and thus the de- 
gree of their asymmetry. General expression for energy 
shift of quantum states of isotope with mass M is sum of 
the normal (virial), specific (the momentum correlation of 
nucleons) and field (nuclear volume) shifts. Symbolically, 
it can be expressed in the form 



AE 



M 



M 



E + E S + E A . 



(8) 



The first term here is the normal shift, being, as seen 
from the expression, large for light element isotopes and 
small for the heavy ones. The specific shift gives essen- 
tial contribution only for the medium weights atoms with 
Z = 10, . . . , 20, and its computation is rather complicated. 
The third term is the nuclear volume shift due to symmet- 
ric charge distribution in the nucleus, depending on the 
integrated product of the mean charge density in nucleus 



and on the square radius of nucleus, i.e. approximately 
holds 

E A oc {p P r 2 M ) , 

where p e is charge density in the nucleus and Tm is cur- 
rent value of the radius. As well known, only the s-state 
electrons have non-zero charge in the centre of nucleus. 
Therefore, the volume shift has largest contribution for 
s-electrons and it is essentially lower for the other states. 

If nucleus is a sphere with radius R and charge distri- 
bution is uniform in it, then 

R 2 = ■ 

For atoms with mass number A > 16 holds approximately 
(in fm units) R = 1.1154A 1 / 3 and thus 



E A cx p e A 2 / 3 . 



(9) 



Naturally, this formula gives only a rough approximation 
for the volume shifts. As seen from the Eq.®, the volume 
shift grows monotonously with atomic weight contrary to 
the normal shift, which diminishes for the heavier nuclei. 

In the isotopes with odd number of nucleons presence 
of magnetic dipole momentum of nucleus and anisotropy 
of the charge distribution generate the atomic energy level 
shifts due to total magnetic dipole and electric quadrupole 
momenta of nucleus, correspondingly. The spectral line 
shifts of the bound states which appear due to magnetic 
dipole, are specified by 



W H 2 = 

and due to quadrupole momentum 

+ 1) -/(/ + !)- J(J- 



W. 



EA 



B- 



21(21 - 1)2J(2J- 1) 



where 



C = F(F + 1) - 1(1 + 1) - J(J + 1) 



(10) 



(11) 



(12) 



In these formulae / and J are the nuclear and total elec- 
tron angular momentum quantum numbers, respectively, 
and F is their vector sum. The quantities A and B de- 
pend correspondingly on the magnetic dipole and electric 
quadrupole momenta of the atomic nucleus. 

4. Separation of mercury isotopes in CP stellar at- 
mospheres 

Generation of elemental and isotopic peculiarities in 
the atmospheres of CP stars is described by the following 
scenario. During stellar evolution line-rich heavy metals 
are pushed to the atmosphere from stellar interiors due to 
the large expelling ra diative force which can exc eed the 
gravity by several dex (jMichaud and Richer! . 120081 ) . How- 
ever, this radiative drive leaves solar isotope ratios unal- 
tered. 
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Isotope separation takes place in stellar atmospheres, 
where the local diffusive radiative transfer is replaced by 
non-local radiative transfer. This favours generation of 
asymmetrical radiative flux in profiles of spectral lines. 
Flux asymmetry within spectral line profiles of atoms and 
ions generates light-induced drift, which is not random 
but cumulative for isotopes, because isotopic splitting of 
spectral lines creates systematic asymmetry in overlapping 
line profiles. 

Study of evolutionary separation of heavy metals, pre- 
dominantly of Hg, and their isotopes due to LID in the 
atmospheres of CP stars has been an essential topic of our 
investigat ions during about las t decade. Computer code 
SMART dSapar and Poolamad . 120031 ; ISapar et all 120071 ; 



Aret et all |2008), composed primarily for modelling of 



stellar atmospheres and stellar spectra, has been supple- 



mented with additional software blocks for computation of 
evolutionary scenarios of diffusive separation of isotopes of 
chemical elements due to radiative acceleration, LID and 
gravity. Mercury has been chosen for the modelling mainly 
because of two reasons. First, it is a typical isotope-rich 
heavy metal. Second, mercury is highly overabundant in 
HgMn stars and its isotopic mixture has been determined 
from observations for more than 30 stars f Woolf and Lam- 
bert. Il999l ; iDolk et al.l . l2003l ). 

It has turned out that formation of anomalous isotope 
abundances in stellar atmospheres, including the domi- 
nance of the heaviest isotope, cannot be explained by the 
diffusion theory without including LID mechanism. 

Since LID is sensitive to the shape of spectral lines, syn- 
thetic spectra must be computed with high resolution at 
all layers of model atmosphere. Radiative flux in optically 




Figure 1: Evolutionary changes of mercury isotope concentrations relative to the initial ones C; = Pi/p® in model atmosphere with T e g 
10 750 K, log g = 4, V TO t = 0, Vt UI b = and initial Hg abundance solar + 5 dex with solar system isotope ratios. 
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Figure 2: Equilibrium mercury concentrations d = pi /pi in model 
atmosphere with microturbulence. Microturbulcnt diffusion coef- 
ficient is assumed to be 50 times the atomic diffusion coefficient. 
Model parameters: T e gp = 10 750 K, logg = 4, p° = solar + 5 dex. 



thin medium of stellar atmosphere has to be found resolv- 
ing the radiative transfer equation in detail. This is more 
difficult and time-consuming than computing the flux ac- 
cording to the so-called diffusion approximation, valid in 
optically thick media of stellar interiors. Precise data on 
hyperfine and isotopic splitting of spectral lines and cross- 
sections of various physical processes are also necessary in 
LID computations. 

Spectral line data for mercury have been compiled us- 
ing different sources and improved by adding isotopic split- 
ting to all available Hg lines . Most of the line data have 
been taken from the lists bv lKurucz and Bell (1 19951) an d 



Vienna atomic line database VALD (jKupka et all 2000l ) . 
Hyperfine and isotopic structure of H g lines and their oscil - 
lato r strengths av ailable in papers bv lProffitt etaD (|l999l ) 



and ISmith (1997) were used to improve our line list. For 
other lines isot opic splitting was calcula t ed us ing relative 
shifts found by IStriganov and Dontsov ( 1955 ) , scaled to 
units [202 — 200] = 1 , giving for other splitted line com- 
ponents [198 - 200] = -0.94, [199 - 200] = -0.80, [201 - 
200] =0.30, [204 - 200] = 1.98. For scaling these relative 
shifts to wavenumbers, we have taken into account that 
[202 - 200] = 0.179 cm" 1 for Hg I, [202 - 200] = 0.508 cm- 1 
for Hg II and [202-200] = 0.600 cm" 1 for Hg III. Compiled 
line list contains about 700 resonance and low excitation 
spectral lines in the wavelength range approximately 800- 
12,000 A for Hg I, Hg II and Hg III, i. e. for ion species, 
which are most important for LID. 

In this paper we present some results of modelling the 
evolutionary separation of mercury isotopes and discuss 
some problems on the way of further studies. The prob- 
lems of LID modelling in the quiescent atmospheres of 
slowly rotating CP stars involve the computation of line 
strengths, formation of line profiles due to different micro- 
physical interaction processes, but also formation of mi- 
croturbulence and weak stellar winds as the phenomena, 
reducing the diffusion. In addition, presence of entangled 
magnetic field can affect the formation of anomalous iso- 
tope abundances, mainly due to Zeeman splitting of spec- 
tral lines. 

Example of a computed evolutionary scenario of sep- 



aration of mercury isotopes due to LID in mercury-rich 
(solar + 5 dex) main sequence stellar atmosphere (T e ff = 
10 750 K, log 5 = 4, F rot = 0, V tmh = 0) is given in Fig.[U 
Macromotions that can suppress diffusion effects have been 
ignored. Consequent frames show relative changes of mer- 
cury isotopic abundances Ci—pij p\ with time throughout 
the atmosphere (—6 < logr < 2). Lighter even-A isotopes 
sink first. Odd- A isotopes are supported longer in the at- 
mosphere due to hyperfine splitting of their spectral lines. 
As stated above, the hyperfine splitting of atomic states 
is essential only for odd-A isotopes. Hyperfine splitting 
mixes the order of isotope lines and makes the picture of 
evolutionary isotope separation more complicated. How- 
ever, as time passes, all lighter isotopes sink and only the 
heaviest isotope 204 Hg remains in the atmosphere. Separa- 
tion in upper layers of the atmosphere proceeds essentially 
faster than in deep layers. We would like to emphasize that 
present evolutionary scenario has been computed for per- 
fectly quiet stellar plasma, without microturbulence, stel- 
lar wind or any other mixing process. Presence of macro- 
motions in the atmosphere drastically slows down the dif- 
fusion and also reduces resulting abundance gradients. 

Final equilibrium isotope separation profiles through- 
out the atmosphere in the presence of microturbulence 
(Fig- HI) have been iteratively computed for the same model 
atmosphere as the evolutionary scenario presented in Fig.[T] 
Microturbulent diffusion coefficient has been assumed to 
exceed the atomic diffusion coefficient 50 times. The heavi- 
est isotope Hg is also strongly dominant, however abun- 
dance gradients are not as steep as in Fig. [1] 

Light-induced drift is generated when the radiative flux 
is asymmetric within the line widths of the light-absorbing 
ion. This asymmetry mainly appears due to line blendings, 
which are usually random. Thus, a direction and value of 
the LID in each spectral line is also random and total 
LID, obtained summing over all spectral lines of studied 
element, is minor. Total value of LID can be large only if 
there is a systematic asymmetry of spectral lines, what is 
the case for the overlapping isotopic lines. Thereby LID of 
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Figure 3: Effect of line blendings on LID of Hg isotopes. Accelera- 
tions of the Hg isotopes due to LID at the first time step (i.e. with 
homogeneous Hg abundance) are shown in the modified logarithmic 
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+ 1 J . Thick lines: only Hg iso- 



topic lines, no blendings with lines of other elements; thin lines: lines 
of all elements are present. 
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mercury isotopes is governed by mutual overlaps of their 
spectral lines, while blendings with lines of other elements 
are less important, because they statistically cancel out. 
Calculations of LID due to overlapping lines of mercury 
isotopes in the presence and without lines of other elements 
confirm this conclusion (Fig. [3]). Blendings with lines of 
other elements certainly affect values of a LID , but do not 
change the main regularities defined by mutual overlaps of 
Hg isotopic lines. 
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5. Hydrogenic line profiles 

Although details of mercury lines are crucial for com- 
putation of LID acceleration, also profiles of all strong 
lines, which overlap with mercury lines, should be com- 
puted with high precision to obtain more realistic picture 
of isotope separation. This concerns also the hydroge n 



ydrogen 
., 2006), 



spectral line series. As we have shown (jSapar et al.l 
the convolution of Holtsmark, Lorentz and Doppler pro- 
files in the result of two analytical integrations reduces to 
the sum of three contribution terms, namely 



*(/3)=xG8)+AG9) + A(/?), 



where 



MP) 



A(/3) 



2L 



4D 2 



and 



0x sin(/3x)e(x)dx, 



x cos(/3x)s(x)dx, 



x 2 cos(f3x)e(x)dx 



e(x) = cxp (-x 3/2 - Lx - D 2 x 2 



(13) 



(14) 



(15) 



(16) 



(17) 



In these formulae j3 = {y — vo)/Ai/s, L — Av^/ Avs and 
D = /S.vn /2/S.vs- Here Avjj, Avl and Av$ are the widths 
parameters of Doppler, Lorentz and Holtsmark line pro- 
files, respectively. 

Another possibility is numerical integration of the con- 
volution of analytical high-precision ap proximation expres - 
sions of Holtsmark and Voigt functions (jSapar et al.Ll200a ). 
The asymptotic series expansions for the small and large 
values of argument (3 have been found by us as a contri- 
bution to this approach. For the intermediate region the 
exponential term can be expanded into series followed by 
integration over circle. 

For main sequence stars Holtsmark distribution shielded 
by the Debye radius, i.e. the Ecker correction is not es- 
sential, because the value of screening parameter, approx- 
imately estimated by 



where Rd is the Debye radius 



/ kT \ 1/2 



(18) 



(19) 



Figure 4: Number of charged particles in the Debye sphere. Large 
Ecker coefficient values S > 50 show that Holtsmark distribution is 
valid throughout stellar atmospheres. 



and N is the number density of charged particles, has large 
values up to deep atmospheric layers. This is illustrated 
in Fig. [U As known, Holtsmark distribution holds well if 
6, being the number of particles inside the Debye sphere, 
exceeds about 50. 

As mentioned earlier, an important problem for mod- 
elling of stellar atmospheres is to find high-precision profile 
function for hydrogen lines. T he Model Microfield Method 
(M MM) initially developed by Frisch and Brissaudl (|l97lh 
and iBrissaud and Frisch dl971 ) a nd further elaborated by 
Stehle et al. (e.g. IStehleL Il990l : IStehle and Jacquemotj . 
19931 ) takes statistically into account also the ion dynamics 
effects. Tables for th e Lym an and Bal mer (IStehlel Il994l ; 
Stehle and Hutcheonl. Il999h. Paschen dStehlel Il9961) and 
Brackett lines ( Stehle and Fouquetl . 2009 ) have been cal- 
culated and are available for astrophysical applications. 
The main difference between MMM and Holtsmark profiles 
lies in the line centre. The profile function obtained with 
MMM has a non-zero plateau at small argument values 
while the Holtsmark distribution approaches zero. How- 
ever, the total contribution of blendings with hydrogen 
lines to the LID of mercury isotopes turned out to be mi- 
nor and eventually also the MMM correction does not es- 
sentially affect our results. 

Splitting of H spectral lines due to the linear Stark 
effect, where the number of splitted line components grows 
rapidly towards the higher members of the series, also must 
be further taken into account. This problem is connected 
with determination of dissolution rates of high Rydberg 
lines. The state survival fractions have been specified at 
least in three different ways, giving the expressions in the 
form 



W n = CN k 



(20) 



where N is the number density of charged particles and 
n is the current value of the principal quantum number. 
These methods are: 



Inglis-Teller overlap cutoff (W n = C±N x n 



the nearest neighbour dissolution (W n = Ci N 1 n 6 ) 
or exponential cutoff, and 

the red-wing Stark components dissolution 



Different partition functions and their dependence on TV 
correspond to these cutoff methods. 

We have used the last expression not only for H line 
series, but for smooth transition from Rydberg line series 
to corresponding continua of all elements. Thus we ob- 
tain smooth transition from line series to the correspond- 
ing continua, but the state survival fractions are probably 
somewhat overestimated. It is extremely difficult to de- 
scribe the real situation, because also a complicated pic- 
ture of Stark splitting of spectral lines must be taken into 
account. Although these problems are weakly connected 
with the LID computations, they are important for mod- 
elling of stellar atmospheres and spectra. 

6. Concluding remarks 

Light-induced drift enables to explain and to model 
evolutionary generation of drastic isotopic anomalies in 
quiescent CP stellar atmospheres, i.e. if stellar wind, merid- 
ional circulation and macroturbulence are lacking. How- 
ever, microturbulence, which is always present in the at- 
mospheres, diminishes essentially the equilibrium gradi- 
ents of isotope ratios. Nevertheless, it has been succeeded 
to explain the main regularity - formation of essentially 
higher abundances of the heavier isotopes of heavy metals 
on the example of mercury. 

Further refinements of a number of physical process 
characteristics are needed. The line profile functions and 
shifts, state survival fractions, impact cross-sections, Stark 
and hyperfine splitting of spectral lines and several other 
quantities should be improved to model evolutionary sce- 
narios of isotope separation in CP stellar atmospheres. 
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